Novel Bound States Treatment of the Two Dimensional 
Schrodinger Equation with Pseudocentral Potentials Plus 
Multiparameter Noncentral Potential 

Metin Akta§* 

Department of Physics, Faculty of Science 
Karatekin University of Qankiri, 18100, Qankiri, Turkey 

September 6, 2011 



Abstract 

By converting the rectangular basis potential V(x,y) into the form as V(r) + V(r,(p) 
described by the pseudo central plus noncentral potential, particular solutions of the 
two dimensional Schrodinger equation in plane-polar coordinates have been carried out 
through the analytic approaching technique of the Nikiforov and Uvarov (NUT). Both 
the exact bound state energy spectra and the corresponding bound state wavefunctions 
of the complete system are determined explicitly and in closed forms. Our presented 
results are identical to those of the previous works and they may also be useful for in- 
vestigation and analysis of structural characteristics in a variety of quantum systems. 
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1 Introduction 



The pioneering proposal of Hartmann [1] namely as the ring-shaped potential and its appli- 
cations [2] has offered much attention in recent times. The so-called Hartmann ring-shaped 
potential V(r, tp) exhibits a special case for the noncentral class potentials. It can be expressed 
as the combinations of an attractive Coulombic term (A/r) and a repulsive noncentral term 
(B/r 2 sin 2 (p) together with the parameters A = 2ao£o^c 2 and B = —qaoEorfa 2 . Here also 
Oo = h 2 /me 2 , £o = —me 4 /2h 2 stand for the Bohr radius, the ground states energy of the 
hydrogen atom respectively with the dimensionless parameters q, r\ and a varying from 1 up 
to 10. 

By proving the evidence of 'accidental' degeneracies in quantum levels, the dynamical invari- 
ance of algebra for the Hartmann potential was firstly handled by the authors Kibler, Negadi 
[3], Gerry [4] and Kibler, Winternitz [5] respectively. On behalf of the work [5], they have 
realized that the dynamical symmetry of quantum systems gives rise to the 'accidental' de- 
generacy which is exhibited the discrete spectrum. Zhedanov [6] also proposed the dynamical 
group symmetry approach to prove the 'hidden' symmetry algebra of the Hartmann and os- 
cillator types ring-shaped potential systems. 

Two distinct forms of noncentral class potentials are realized in the literature. The first form 
refers to the Hartmann type ring-shaped potential. The second form firstly introduced by 
Quesne [7] to investigate the 'accidental' degeneracy corresponds to the oscillator type ring- 
shaped potential. It is employed by combining the terms of harmonic oscillator (~ r 2 ) and 
(r, y?)-dependent noncentral potential [8]. It is pointed out that in the first classes the lead- 
ing term implies the Coulombic discrete energies though in the second ones the leading term 
suggests the harmonic oscillator basis discrete energies. 

Quite a few applications have been encountered in such works on the subject of noncentral 
potentials (NCPs) from past to the present time. For example, general relations between the 
elastic constants and the central forces in hexagonal materials as well as noncentral forces 
in fee monatomic structures were given by Johnson [9]. Some certain structures with these 
type potentials were also studied by the authors [10-16]. Ermakov type invariants were also 
discussed with respect to them by Makowski in exact classical limit of quantum mechanics 
[17]. Besides, the derivations were achieved for determining the stress and elastic constants 
in systems of particles via noncentral two-body potentials by Murat and Kantor [18]. Several 
authors has employed them for the aspect of scattering analysis in continuum bound states 
[19-25]. Fredholm theory was applied to the Lippmann-Schwinger equation and the general- 
ized Levinson theorem was also proved for NCPs [19]. Forward scattering in a system with 
those potentials was discussed by [20]. 

Recently, much considerable effort for variety forms of NCPs has been expanded on the solu- 
tions of Schrodinger, Dirac and Klein-Gordon equations. The Feynman's path integral treat- 
ment [26-30] and the Green's function technique [31, 32], the (Lie) algebraic/group theoret- 
ical approach [33-36], nonbijective canonical transformation [3, 37], supersymmetric (SUSY) 
quantum mechanical formalism [38-46] and the NU-analytic method [47-58] as well as the 
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applications for both relativistic [59-71] and other nonrelativistic [72-82] cases are available in 
the literature. 

Due to the above reasons, NCPs may be helpful for determining the structural properties 
of systems in such cases (e.g. elasticity and stress factors, point-defects of surfaces etc.) in 
variety fields of physics and chemistry. The ring-shaped like structures, i.e., cyclic polyenes, 
benzene and benzene-like structures (e.g. graphene), and interactions between deformed pair 
of nuclei are good examples for the discussion. 

The key idea 'noncentrality concept' of all potentials plays a privilege role to provide the 
extensive solutions for concerning quantum systems. That is, if one of NCPs preserves the 
separability condition for achieving the solutions of system, it can therefore be split into two 
parts for such cases of radial and angular dependent equations that they may readily be solved. 

Our primary goal of this work is to present the expressions of the exact bound states both 
the energy spectra and the wavefunctions of the two dimensional Schrodinger equation by 
proposing two definite potential forms in which they are called the pseudocentral potentials 
plus multiparameter NCP. It is the author's hope to further contribute to applying our ap- 
proach and results for the most recent quantum systems and for exactly solvable systems 
with centrally style potentials [83, 84]. It may also be applicable for conforming the initial 
boundary- value problems of certain physical systems (e.g. Dirichlet, Neumann etc.) in two- 
dimension [85-87]. 

The organization scheme of the study is given as follows: Section 2 is devoted to construction 
of the two different rectangular potentials and converting of them to be solvable form for the 
Schrodinger equation in plane polar coordinates. Section 3 covers the solutions of separated 
equations via the NU-method. Section 4 is responsible for dealing with the conclusions and 
remarks. 

2 Potential Cases and Separation of Variables 

Here we will firstly describe the essential steps of converting the rectangular basis potentials 
as Vj(x : y) and Vjj(x,y) into the plane polar forms as V/(r) + Vj(r : ip) and V n (r) + V n (r,(p). 
Secondly, by applying the separational procedure to the two dimensional Schrodinger equation 
(2-dim SE) then we are going to deal with the solutions of separated equations. 

We propose the rectangular form physical potentials as follows 

--Up**)'-* 



and 



Bx 4 + cy 

(x 2 + y 2 ) 



+ Dx 2 + Fy 2 + G(x 2 + y 2 ) 



(1) 
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V n (x,y) = (A iy fx 2 + y 2 + 



B l 



\Jx 2 + y 2 J x 2 y 2 



Bx 4 + Cy 4 
(x 2 + y 2 ) 



+ Dx 2 + Fy 2 + G(x 2 + y 2 



(2) 



The squared leading terms in equations (1) and (2) refer to the pseudo Coulombic/modified 
Kratzer and the pseudoharmonic circular oscillator potentials as well. The rest parts of them 
correspond to the equivalent multiparameter potentials which will be converted into the non- 
central case. 

Let us consider the cartesian coordinate transformations as [88] x = rcosy?, y = rsimp, 
(p = tan~ 1 (y /x) and r = ^/ x 2 + y 2 , by putting these into the above equations then they are 
readily converted into the plane polar forms as 

V(r) = V I (i) + V I (T,(p) = (^ + ^ + ?f) + ^{BcotV + CtanV + DcscV 

+ F sec 2 (f + G sec 2 ip esc 2 </?j (3) 



with the constants kq = B 2 = D e , K\ = 2A B = 2D e r e and k^ = A\ = D e r 2 as in Ref. [50]. 
The consecutive terms, in parenthesis, refer to the pseudocentral potential (PCP1) called as 
the pseudo Coulombic/modified Kratzer case and the rest part constitutes the multiparameter 
noncentral potential labelled as (mp — NCP). By applying the conversion process just only 
the equation (2) gives 



V(t) = V n (r) + V n (v, <p) = V + A{ r 2 + 



Bl 



(mp - NCP) 



(4) 



where we have used the notations V = 2A 1 B 1 = (k/A)tq with A 1 = (k/8) 1 ^ 2 , and 
Bi = (k/8) 1 / 2 t-q as [84]. The consecutive three-term, in equation (4), establishes the pseu- 
docentral harmonic oscillator potential (PCP2). 

Let us now regard the two dimensional rectangular basis Schrodinger equation (SE). It can 
be written 



<9 2 ^ d 2 ^ 2m 



dx 2 



(5) 



By employing the task of the polar coordinate transformation on equation (5), we rewrite it as 
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Id ( d\ Id 2 2m rrn . ,.' 



tt(r,<p) = 0. 



(6) 



The separation of variables procedure cannot be carried out when rectangular coordinates are 
employed since the potential energies (l)and (2) are functions of (x 2 + y 2 )~ 1 and (x 2 + y 2 )~ 1 ^ 2 ■ 
Due to the fact that they are not split into terms, it is required changing to the plane polar 
coordinates. 



Using the ansatz wavefunction in equation (6) for successively concerning the equations (3) 
and (4), 



M>(r,y9) = ^U(r)$(y?) 



(7) 



leads to the following equations in which two of them are radial-dependent correspond to the 
first and second classes of NCPs 



(8) 



and 



d 2 U E 
dr 2 



+ E 



At 2 --) U b (t) = 0, 



(9) 



and the angle-dependent equation 



d 2 $ 
d^ 2 



M 2 - 



'D + Bcos 2 ip\ /F + Csin 2 v?\ 



sin 2 ip 



cos 2 tp 



a 

sin 2 (p cos 2 cp 



${<p)=0. (10) 



In these equations, we use the specifications E H rs = 2m(E — kq)/Ti 2 , A = {2mKi)/Ti 2 
and A = [R 2 + (M 2 - \)] with R 2 = (2mK 2 )/h 2 = (2mA 2 ) /h 2 ; E RSO = 2m(E - V )/h 2 , 
A = [2mA 2 ) /h 2 , and r = + (M 2 - |)] with B 1 = (2mB 2 /h 2 ). In addition, E HRS and 
Ei?5o belong to the energies of the Hartmann ring-shaped (HRS) and the ring-shaped oscil- 
lator (RSO) potentials, respectively. 



Besides we label the energy parameters for noncentral potential part as B = (2mB)/h 2 , 
C = (2mC)/h 2 , D = (2mD)/h 2 , F = (2mF)/h 2 and G = {2mG)/h 2 with the separation 
constant M recognized as the angular momentum. 
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3 Achieving the Solutions of Separated Equations 
3.1 Solution of the <p— Angle Dependent Equation 

Let us begin firstly to examine some solutions of the equation (10). By introducing a trans- 
formation as t = sin 2 ip hence the transformed hypergeometric equation holds a form as 

d 2 $ (l-2t)rf$ 1 f 2 n W/N « 



where the parameters are a = (-M 2 + B - C), (3 = (M 2 — D — F) and 7 = -(£ + D + G). 
By following [47], a hypergeometric equation is described by 

d 2 u fit) du ait) . . . . 

"TT + -7T^ + — TTX u {t) = 0, 12 



then comparing the equations (11) and (12) term-by-term allow us to write the polynomials 
f(t) = (l-2t) <r(i) = 2t(l - and cr(i) = (at 2 + + 7). (13) 

They are all used in solving of the quadratic equation 

TT(f) = ^ (a' - f) ± i^-f^+^-flr), (14) 



where the prime factor of a denotes the differential at first degree. In this equation the deter- 
mination of k is very essential step for the calculation of ir(t) in which the inner part of the 
square root is required to become the square form with respect to the polynomials. 

The substitution process of polynomials for n(t) and the arrangement of it for k provide us 



7r(t) = -(l-2t)±-yj ai t 2 + a 2 t + a 3 (15) 



with the parameters a± = (—4a — 8k + 4), a 2 = (—4/3 + 8A; — 4) and as = (1 — 47). The 
values of k which bring about the releasing from the square root of n(t) can be determined 
by solving the quadratic equation 
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k 2 + (-p - 2 7 - -) k + ^ (f3 2 + 2(3 + a + 4 7 (1 - a)) = 0. 



(16) 



Then it yields the double roots as follows 



fei,2 = (t + § + \) ± ^ -47) {1-4 [ 7 + (« + /?)]}• 



:i7) 



By inserting these roots into the equation (15) leads to the fourfold roots for n(t) as 



tt 1)2 (0 = 1(1 - 2t) ± \ { [yi - 4 7 - 71 - 4 [7 + (a + /?)]] t - ^/l - 4 7 } (18) 



for negative root of k — k\ and 



7r 3>4 (i) = \{l - 2t) ± \ { [v/l - 4 7 + VI - 4 [7 + (a + /?)]] t - ^1 - 4 7 } 



(19) 



for positive root oik = k2, respectively. The resulting equations (13), (17), (18) and (19) play 
a key role to establish both the energy spectra and the wavefunctions of the system. 

3.2 Bound States Eigenvalues 

Now we can put forward the procedure for determining the energy spectra of the systems. For 
the purpose of this, we again follow [47] and a novel form of the energy equation is adapted 
as 



A = A„ = k + 



dn 
~dt 

( d 



-n I — (r + 2vr) - 2n + 2 



(20) 



Each root k± for ^^(i) and k2 for (t) can be applied which one of them is satisfied the 
condition r' < \t' = J^(f + 2ir)]. Possible solutions have often been carried out by the 
negative roots of the functions n(t). Hence we accomplish the result of the angular momentum 
M which implies the reasonable meaning for physically 
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M 



[D + F + 2(B + G)} + - -1 + y/[l + A(B + D + G)} [1 + 4(C + F + G)] 



+ |(2n + l) (2n + 1) + -(y/l + A(B + D + G) + ^1 + A(C + F + G)) + lj(21) 
with n = 0, 1, 2 . . . 

One of our basic interests in this step is to show and present how to get the energy eigenvalue 
results for the radial-dependent equations (8) and (9) with the potential functions (3) and (4) 
respectively by employing the analogy procedure. It is obvious that the equation (38) and its 
converted equation (50) given as [48] are identical form to the equation (8). Accordingly, it 
should also admit to be similar form solutions. Following the same procedure as in [48] and 
accepting the energy parameters of [50] give rise to 



Eni — «o 



8m 



ni + ^j +^Vl + 4A 



2m 9 9 
n T)r 

e ^2 e 1 e 




2mD e rl ~ 
c + M 



-2 



m = o,i,... (22) 



Note that the equation (57) in [48] and the energy expression (22) referred to (PCP1) plus 
(rap — NCP) potential exhibit similar characteristic by virtue of the presence of the Coulom- 
bic term in potential energy. 



The energy spectra expression for (PCP2) plus (rap — NCP) can also be determined anal- 
ogously by dealing with [56] and [84]. The equation (15) and its converted form (18) in [56] 
hold similar structure to the form of the equation given by (9). Following the procedure as 
[56] step-by-step and the arrangement of the result allow us to write analogously 



E 



"2 



Vn — Twj 



1 



(2n 2 + l) + -v / l + 4r 



i 2 h [k 

-nr - -J — 



(2n 2 + 1 




n 2 = 0, 1, 



(23) 



where k = Arau 2 with to = y (2raA\)/h 2 and the angular momentum M. 
3.3 Bound State Wavefunctions 

The essential steps for performing a set of wavefunction for the complete system should now 
be examined. In this point we are going to introduce the proper notations in the solutions 
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of/c = A;2 = 7+ f + i (\/ (1 + (1 + (3 2 ) ) for which this is the correspondence polyno- 
mial of ir(t) = tt 4 (0 = |(1 - 2t) - \ { [v/TTA + v / TT^] t - y/TTWi} with ft = -4 7 and 
/3 2 = — 4 [7 + (a + /?)] as well as a, /3 and 7 are given in the equation (11). 

Let us start to write the equation 

and from the equation (13), straightforward calculations give for 

*(t) = [t(l - i)] 5 / 2 (25) 

where 5 = [1 + (VI + A + VI + #2 )/4]. The weight function [47] g(t) is determined by 

j t (ga) = gr (26) 

with the expression r = f + 2?t. The calculating procedure yields 

= [*(1 " t)} 6 - 1 (27) 

where 5 — 1 = (VI + A + VI + ^2 )/4. Then the Rodriguez formula is 

*<«> = ^SF^'M- (28) 

where C n is the normalization constant. Use of p(t) and a it) provides 

y n (t) = C n P^^(t), t = sinV (29) 



with C n = 2 n C n , /ii = VI + Ji and fi 2 = VI + fa as well as Pjf* 1 ' ^{t) stands for the Jacobi 
polynomial. The net wave function can thus be written as 



$(t) = <f>(t)y n (t) 

= C no [t(l-t)f 2 P^^(t) (30) 



In order to construct the wavef unctions for the radial-dependent equations (8) and (9), we 
should consider the equations from (58) to (62) as [48] and the equations from (35) to (42) as 
[56], respectively. By means of the analogy procedure one gets 

U(t) = Ua(t) — 0(r) y n ( v ) 

= ^r^e-^^^r) (31) 

and 

U(s) = U B (s) = 0(s) y n (s) 

= C n2 s^e-^ s K 2 (s), s = r 2 (32) 

where we have used the short notations p — 1 + 2y/l + 4A, ( = 1 + y/1 + 4T and q = (£ — 1)/2. 
Also (r) and L^ 2 (s) stand for the Laguerre polynomials, respectively. 

3.4 Complete Solutions of Bound Wavefunctions 

In this part it is required to gather all wavefunctions in which they are satisfied by the equa- 
tions of (8), (9) and (10) with respect to the ansatz wavefunction (7). If we substitute the 
equations (30) and (31) into the equation (7) one yields 

tf(r,p) = N r (p ~5) e -VE M s r ( sin2 ^) 2<5 p(£i. ^)(sin 2 </?) (33) 

where N = (C no ■ C ni )/2 S , sin2y9 = 2 sin ip cos ip and no, n\ = 0,1,2.... By employing the 
cartesian coordinate transformations then it is transformed into 

V(x, y) = N (x 2 + y 2 )^)/ 2 e -VE^~s [ xy /^ + y 2 )] 25 

In this time we consider the equations (30) and (32) for ansatz equation (7) one gets 

tf(r,p) = N t^/ 2 r 2 L « 2 (r 2 ) (sin2y9) 25 P^ 1 ' " 2) (sin 2 y9) (35) 

with n , n 2 — 0, 1, 2 Consequently, it can be converted into 
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*(x,y) = N (x' + y 2 )^ 4 e- 1 ^^ 2 ^ [xy/(x 2 + y 2 )] 25 

x K 2 (x 2 + y 2 )P { :;^ 2) [y 2 /(x 2 + y% (36) 

where N = (C m ■ C n2 )/2 S . 

4 Conclusions and Remarks 

In this study we have performed the exact analytical bound state solutions both the energy 
spectra and the corresponding wavefunctions for the two dimensional Schrodinger equation in 
plane polar coordinates. The rectangular basis potentials proposed by the equations (1) and 
(2) are converted into the planar forms in solving of the 2-dim SE through the NU-analytic 
technique. 

Let us remark that the discussion of noncentral class potentials for such quantum systems is 
limited in the literature. This is the main reason of all not satisfying the condition of ana- 
lytical solvability. On the other hand, the radial-angular dependency of such potentials will 
provide us for examining the structural properties of planar quantum systems, i. e. elasticity, 
stress factor and point-defects of surfaces etc. It is expected that our straightforward approach 
may therefore offers a solution or provide a model for other identical systems described by [84]. 

With regard to a variety of potential parameters, possible cases are briefly examined for en- 
ergy expressions presented by (22) and (23). Certain examples are given as follows: 

Letting B = C = D = F = G = the angular-dependent equation (10) and the equation (21) 
are ignored, therefore this just refers to the energies of the pseudo Coulombic potential (PCP1) 
and the pseudoharmonic oscillator potential (PCP2). For D^O, and B — C — F — G — 0, 
we accomplish the energy eigenvalue results of the pseudo Coulombic type Hartmann ring- 
shaped potential and the pseudoharmonic oscillator type ring-shaped potential as well. 

As a final remark we conclude that the two dimensional angular momentum operator defined 
as [88] 

m> = -%h— 

dip 

with the eigenvalue (KM) 

= (hM)^ 

where the \I/ corresponds to the complete solutions of bound wavefunctions (33) and (35), 
angular momentum expression of our proposal system is also defined by (21). 
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